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CN Abstract 

^ A new and rather broad class of stationary (i.e. stochastically translation in- 

•"^ variant) random tessellations of the d-dimensional Euclidean space is introduced, 

which are called shape-driven nested Markov tessellations. Locally, these tessella- 
tions are constructed by means of a spatio-temporal random recursive split dynamics 
governed by a family of Markovian split kernel, generalizing thereby the - by now 
classical - construction of iteration stable random tessellations. By providing an 
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sumptions on the split kernels (shape-driven), there exists a unique time-consistent 
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the typical cell and some aspects of the first-order geometry of these tessellations are 

^ in the focus of our discussion. 
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1 Introduction 

Modern stochastic geometry has a rapidly growing demand for non-trivial and flexible, yet 
mathematically tractable models for spatial random structures of the d-dimensional Eu- 
clidean space R*^, d> 2. The need for such new model arises in particular in the theory of 
random tessellations, as they are widely used for modeling structures ranging from classical 
geology to the more recent telecommunication science (cf. [17, 23]). Stationary random 
tessellations that are stable under the operation of iteration - so-called STIT tessellations 
- were introduced recently by Mecke, Nagel and WeiB [11, 12, 14] and they may serve as 
a new mathematical reference model beside the classical Poisson hyperplane or Poisson- 
Voronoi tessellations. On the other hand, the first-order geometry of these tessellations is 
already determined by two parameters, the surface density and the so-called directional 
distribution. Moreover, STIT tessellations have Poisson-typical cells, which is to say that 
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the distribution of the typical cell of a STIT tessellations coincides with that of a Poisson 
hyperplane tessellation having the same parameters. Thus, the size (measured for example 
by some of the intrinsic volumes) of the cells of STIT tessellations is rather inhomogeneous 
and varies from arbitrarily small to rather large. However, tessellations with such behavior 
are often not suitable for practical applications. For these reasons, it is desirable to intro- 
duce more flexibility into the model. One could think for example of a hard-core condition, 
which excludes in some sense 'small' cells or of a condition that keeps the size of large cells 
under control. In the Poisson- Voronoi reps, the Poisson-Delaunay case a related theory 
was developed in [4, 5]. A first attempt in the case of iteration stable and related models 
was recently made by Cowan [2], where a number of random cell splitting mechanisms of 
finite volumes are discussed, which are mostly independent of the cell geometry (so-called 
geometry independent apportionments). However, neither the problem of existence and 
uniqueness of a related whole-space tessellation nor main properties of such tessellations 
(such as the typical cell distribution or scaling relations for example) have been established 
in that paper and remained open questions. 

In the present work, we start by describing a general class of what we call nested tessel- 
lations. Locally, their construction can be interpreted as outcome of a random process 
of cell division, which may roughly be described as follows: Each cell is provided with a 
random life-time, which is related to the geometry of the cell, to its environment and even 
to the running construction time. Once this life-times is expired, a random hyperplane 
hitting the cell is chosen according to some distribution, which is allowed to depend on the 
cell's geometry, its environment and also on time. It is introduced in the cell, cut-off by 
the cell's boundary and divides it into two polyhedral sub-cells. The whole construction 
is continued independently and recursively in both of the sub-cells, in such a way that 
the respectively new introduced hyperplanes are also cut-off by the boundaries of their 
mother-cells. Formally, the construction is described by a family of Markovian split-kernel 
and uses the general framework of pure jump Markov processes. In [19] it was shown that 
the same informal description can also be provided for the infinite volume construction. 
Thus, each general nested tessellation admits a unique family of split kernels. Naturally, 
it is of particular interest if also the converse is also true, i.e. whether a given family of 
split kernels determines a nested tessellation and, moreover, if this tessellation is unique 
and 'well behaved'. It seems that problem, in its general form is very difficult to handle. 
It is the main purpose of the present paper to prove a result of this kind for the class of 
what we call shape-driven nested Markov tessellations. They are nested tessellations in the 
above sense, whose characterizing split kernels satisfy a number of regularity conditions. 
Beside the existence and uniqueness theorem, we also investigate the geometry all these 
so-constructed tessellations have in common. In the planar case, the most explicit results 
are available, whereas for d > 3, mean values depend on an increasing number of param- 
eters, determined by the precise cell splitting mechanism. For this reason and because 
of the rapidly increasing geometric complexity of the structures, most of our mean value 
formulas will be restricted to the practically relevant cases d — 2 and d — 3. Moreover, we 
discuss in detail the geometry of the typical cell, which eventually leads to a simulation 
method via Monte-Carlo techniques. 
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The paper is structured as follows: In Section 2 we recall in some detail the - by now 
well known- description of STIT tessellations and figure-out the main construction prin- 
ciple underlying the structure of general nested tessellations in terms of their so-called 
spht kernels. After that, in Section 3, the class of shape-driven Markov tessellations is 
introduced by assuming certain, yet flexible, regularity conditions of the split kernels and 
the principal existence and uniqueness result is formulated. Section 4 is focused on general 
mean values for these tessellations in terms of parameters determined by the cell spliting 
mechanism. The argument of our main theorem as well as a global construction of the 
tessellations under consideration is the content of Section 5. 



2 Iteration Stable and General Nested Tessellations 

We consider a translation-invariant measure A on the set % oi {d — l)-dimcnsional (affine) 
hyperplanes in with d > 2 and note that under the standard identification of a hy- 
pcrplanc H = {x E M'^, = r}, u G §d-i, r > with the pair {u,r) G S^-i x ]R+, 

the translational invariance leads to the decomposition K = pTZ ® £+, where is the 
Lebesgue measure on R+, TZ is an even probability measure on the unit sphere ^d-\ and 
p G IR+ is the surface density of A. In the particular case where p = 1 and where TZ is 
the uniform distribution on S^-i, the corresponding measure is isometry invariant and is 
denoted by Aiso in this paper. Often, Aiso is referred to as the invariant hyperplane measure. 

An iteration stable random tessellation 1a (t) with driving measure A and time pa- 
rameter t > is constructed according to the following rules, originally due to Mecke, 
Nagel and Weiss, see [11, 12, 14], and hence referred to as the MNW-construction in 
the sequel. We begin at the time with a compact and convex subset W dW^ (called win- 
dow here) having nonempty interior, regard it as the initial cell for the tessellation under 
construction and assign to it an exponentially distributed random lifetime with parameter 
^{[W]), where 

[W\ -.^{H eH, HnW j^^} 

stands for the set of all hyperplanes hitting W. Upon expiry of its lifetime, the cell W 
dies and splits into two sub-cells and W~ separated by a hyperplane in [W], chosen 
according to the law 

A{-)/A{[W]). (1) 

The resulting new cells and W~ are again assigned independent exponentially dis- 
tributed random hfetimes with respective parameters A([iy"'"]) and A([iy~]) (whence ge- 
ometrically smaller cells live stochastically longer) and the entire construction continues 
recursively, until the deterministic time threshold t > is reached. The cell-separating 
{d, — l)-dimcnsional facets arising in subsequent splits are usually referred to as maxi- 
mal polytopes (which are in this paper always {d — l)-dimensional) or I-segments in 
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the planar case d — 2, a,s assuming shapes similar to the hteral /. We will also speak 
about maximal polygons in the spatial case d = 3. The resulting random tessellation of 
FT C M'' is denoted by Y^it, W) and is usually called the STIT tessellation (for STable 
with respect to ITerations) with driving measure A and time threshold t. By noting that 
for any for compact convex V C W the trace the MNW-construction leaves on V C W 
coincides in law with the MNW-process for Y\{t, V), we find that Y{t, •) enjoys a spatial 

consistency property in that Yj^{t, W) r\V ^Yj^{t,V) , where = stands for equality in dis- 
tribution. Note that by facets, cells and regions of the restricted tessellation Y/^(t, W) fl V 
we understand the non-empty intersections of, respectively, facets, cells and regions of 
^a(^, W) with V. The consistency readily allows us, by utilizing the consistency theorem 
[18, Thm. 2.3.1], to construct the whole-space process Y\{t) := yA(^,IR'^) in such a way 
that YA{t, W) = Yj^{t) n W. Note that, by construction, all cells of Yf,{t, W) and its global 
variant Y\{t) are necessarily convex (and even polyhedral if we assume the initial window 
1^ to be a polytope). For terminological convenience we shall use the name region for a 
general cell of the tessellation. 

To justify the name iteration stable for the tessellation Y\{t) we note that if we insert an 
independent copy of the tessellation Y\{s) fl c into each active cell c of the tessellation 
Yj^{t), the resulting iterated tessellation Y\{t)SY\{s), arising as the local superposition 
of i.i.d. copies of Y\(s) within the cells of Y\(t), coincides in law with Y\{t + s). On the 
other hand, it is easily concluded from the MNW-construction that 

lA(tl) = ^lA(t2), (2) 

which is to say YA(ii) coincides with YA.(t2) re-scaled by the factor ^2/^1- Consequently, 

Yj^{t) = mYj^{tm) = m{Y^{t) ffl . . . ffl Y^{t)) =: mYj^itf"^ (3) 

^ V ' 

m 

for t > and m G N. The relation (3) is what usually goes under the name of stability 
with respect to iterations, whence the name STIT. 

The stability relation with respect to iteration is one of the crucial properties for studying 
the MNW-construction. In fact it turns out [14, Corollary 2] that each translation invariant 
random tessellation of positive and finite surface density arising as a re-scaling of its iter- 
ations does admit a MNW-representation as above. Moreover, these are the only positive 
and finite surface density tessellations arising as re-scaled limits of iterations of stationary 
tessellations as shown by Theorem 3 in [14]. We like to point out that the STIT model 
enjoys truly remarkable properties with a lot of non-trivial, yet explicit, results available 
by now and under active development going on, compare with [19, 20, 21, 22, 24, 25]. 

It is convenient from the formal viewpoint to regard the random STIT tessellation Y\ — 
Y\{t) as the random collection of maximal polytopes marked by their respective birth 
times, that is to say a marked point process in the space of {d — l)-dimensional compact 
polytopes in M*^ carrying marks in [0,t]. With probability one this collection 
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[Empty start] is empty at the time 0, 

[Local finiteness] is locally finite in a sense that a compact subset of M"^ intersects only 
a finite number of maximal polytopes, 

[Iterative binary subdivisions] and enjoys the property that each facet (this is a 
maximal polytope) born at time < s < t is contained in a unique convex cell and 
splits it into two convex sub-cells. 

Each collection of time-marked (ci — l)-dimensional polytopes with these three properties is 
what we call a nested tessellation in this paper. Note, that unbounded cells and regions 
are permitted at this stage. The space of all nested tessellations of a given convex set W C 

R'^, possibly with W = R"^, is denoted by Y[o,t](l^) and we abbreviate Y[o,t] := Y[o,t](]R'^). 
Endowing the space of maximal polytopes with the usual Hausdorff distance (see [18, Chap. 
12.3]) allows us to define a natural vague topology on Y[o,t](W^), see [3, Definition A2.3.1], 
which induces the corresponding Borel cx-field on Y[o,i](iy). Note that upon forgetting 
the time marks we may also interpret a tessellation y e Y[o,t](PF) as closed subset of R'' 
arising as the union of all its maximal polytopes. It is easily verified (compare with [18, 
Thm. 12.2.3]) that this yields a continuous map from Y[o^t]{W) to the space J^{W) of closed 
subsets of W with the usual Fell (hit-or-miss) topology [13, Def. 1.1.1]. Another important 
measurable map is the time restriction map Lg ■ Y[o,t](H^) — )■ Y[o,s](H^), < s <t, removing 
from its argument tessellation all maximal polytopes with birth times exceeding s and 
keeping only those born no later than at the time s. Note, that this, yet measurable, 
mapping is in general not continuous, which may be seen from the remarks on page 567 in 
[18]. By a random nested tessellation in C M'^ we shall obviously mean a random 
element taking values in the space Y[o,t](iy) for some t > 0. For a random nested tessellation 
Y (note that we use y for generic deterministic nested tessellations and Y for generic 
random ones) we write 

Y{s) -.^LsiY), 0<s<t, 

which yields tessellation- valued right-continuous stochastic process (y (s))sg[o,t] on Y[o,t](VF) 
or Y[o.t]. We also write 

Y{s,W) = Y{s)nW 

for s G [O.t] and W C M'^. Clearly, (y(s))sg[o,t] is a (possibly non-homogeneous) Markov 
process, because Y{s) uniquely determines Y{s') for all s' < s. In this paper we will iden- 
tify a nested random tessellion Y e Y[o,t] with the process (F(s))5g[o,t] induced by the time 
restriction map. 

To get more insight into the structure of these processes, we restrict our attention to ran- 
dom tessellations with laws locally absolutely continuous with respect to that of Y\.^^{t), 
that is to say C{Y{t, W)) < C{Yji.^^{t, W)) for all compact convex windows W CR"^ (read 
£(...) as the law of ...). Note in particular that this condition is satisfied for y = iff 
A ^ Aiso. For a fixed compact convex window W <Z R'^ and s G [0,t], by the local abso- 
lute continuity as assumed above, the general theory of pure jump Markov processes [6], 
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[10] guarantees that there exists a split kernel ^^^h that for each y e Y[o,5](VF) and 
c e Cclls(y), $^(-|c, is a finite measure on [c], moreover $^ depends measurably on c, 
y and s and, furthermore, $^ is such that the process {Y{s, W) — Y[s) fl H^)sg[o,t] evolves 
in law according to the following recursive split dynamics 

[Recursive split dynamics] 

• Begin at the time with y (0, W) 0. 

• In time s G [0,t], for each c e Cells(F(s, W)), with intensity 

^YidH\c,Y{s))ds, (4) 
split c with hyperplane H e [c] selected according to the law 

<^>Y{dH\c,Y{s)) 

\^Y{[c]\c,Y{s))\ ^ ' 

into two new cells c^{H) and c~{H) (for definiteness say with c^{H) lying in 
the left half-space determined by iJ), thus adding to Y{s^W) a new maximal 
polytope if n c marked with the birth time s. The resulting new tessellation is 
denoted hy y[c(Z)s H]. 

More formally, it means that Y[s, W) is a (non-homogeneous) Markov process with gen- 
erator given by 

[^7f]{y)= E I [f{y[^^sH])-f{y)]^'Y{dH\c,y\ (6) 

ceCells(2/)'^W 

defined for all bounded measurable / : Y[o,t](H^) M. In particular, one can directly 
conclude from (5) in view of (1) that for the STIT tessellation Ia(^) we have 

^Y^^{dH\c,y)^lHe[c]^{dH), (7) 

whence the kernel does not depend on 1^ for c C 1^ and also not on the time parameter 
s and on the environment y. 

In intuitive terms, all this means that a generic finite volume nested tessellation process 
Y{s^W), s G [0,t], admits a representation where each cell splits with rate given by the 
total mass of its split kernel as depending on the surrounding configuration, whereas the 
probability law - arising upon normalising this kernel - governs the splitting mechanism, 
as determining the cutting hyperplane. Importantly, the same informal description is also 
valid in infinite volume. 

Proposition 1 [Recursive split dynamics] There exist split kernels such that for each 
y £ '^[0,3] o-nd c e Cells(y), $5(-|c, y) is a finite measure on the space [c], moreover $s 
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depends measurably on c and y and, furthermore, $s 'is such that the infinite-volume Markov 
process (y(s))sg[o,t] admits generator 

M{y)= Yl f [f{y[c0sH])-f{y)]^,{dH\c,y) (8) 

ceCells(2/) "^W 

defined for all bounded f which are a{Y{s) : s e [0,t]) -measurable for some bounded convex 
W gW^. Moreover, for He [c], 

^Y(dH\c,y) ^'E[^s(dH\Y(s))\Y(s) nW ^ y]. (9) 

Note that the formal statement (8) simply means that, in analogy to the finite volume 
set-up, also the infinite volume process Y{s), s G [0,t] unfolds in time according to the 
[Recursive split dynamics], where each cell c sphts by a hyperplane H with instan- 
taneous intensity ^s{dH\c,y)ds at time s in environment y. The relation (9) states the 
intuitively obvious fact that for a cell c its instantaneous split rate hj H E [c] coincides 
with the whole space split rate by H, averaged over the environment outside W. 
A proof of Proposition 1 in an extended setting is given in [19], but we also give an argu- 
ment in Appendix A in order to keep the paper self-contained. 

Whereas it is relatively easy to show that, as stated in Proposition 1, each nested tessella- 
tion admits a unique family of characterizing split kernels, the converse question whether 
a given family of kernels determines a tessellation and if this tessellation is unique and well 
behaved, is in general very difficult. In the present paper we shall handle this problem for 
a special class of so-called shape-driven Mcirkov kernels, which are introduced in the 
next section. 

3 Shape-Driven Markov Tessellations 
3.1 The Basic Setting 

Proposition 1 provides a very general structural characterization of random nested tessel- 
lations, at least those locally absolutely continuous with respect to the law of the STIT 
tessellation Y\.^^{t). In particular, it indicates - in agreement with ideas recently developed 
by Cowan [2] - that it is natural to consider the splitting mechanism in iterative nested 
tessellations as consisting of two separate ingredients 

• split intensities, corresponding to the total masses |$s([c]|c, of the split kernels, 

• split geometry, corresponding to the normalized kernels $5(-|c, y)/|$5([c]|c, in our 
setting. 

Whereas most of the work so far was concentrated on the first intensity-related point, 
while usually keeping the geometric ingredient of the kernels just homogeneous, compare 



7 



with the so-called geometry independent apportionment considered in [2] , here our interest 
focuses on the second point. Our theory concentrates on a certain subfamily of admissible 
kernels chosen in such a way that the resulting class of tessellations is rich and flexible 
enough, but yet mathematically tractable. We begin with a natural spatial Markovianity 
postulate for the nested tessellation Y: 

[Markov property] Given any two regions Ai and A2 of Y{s), < s < t, with disjoint 
interiors we have Y{s, Ai) and Y{s, A2) conditionally independent. 

This condition simply means that the behavior of the process within a cell - and recur- 
sively, its sub-cells should they arise - does depend on the cell only through its geometry 
and does not depend on what happens outside of it. In terms of split kernels of Proposition 
1 this is clearly equivalent to having the kernel ^s{dH\c,y) depending only on dH and c 
but not y. That is to say, ^s{dH\c,y) = ^s{dH\c) for all environments y, which even 
better justifies the use of notion of Markovianity in this context. 

Once we assume the so-understood spatially Markovian tessellations, it is important to 
note that if we just change the split intensities of the process while keeping its split ge- 
ometry unchanged, that is to say if we multiply $s('|') by scalar values, then we end up 
with another tessellation process amenable to the original, unmodified one by local time 
re-parametrization over individual regions, yet with the overall tessellation geometry re- 
maining the same upon forgetting the time marks. This motivates us to stick to the usual 
canonical split intensities in this work, which enjoy particularly good properties due to 
their natural relationship to the previously developed STIT theory. Namely, in the sequel 
wc shall always assume that the split intensity is given by the A-mcan width (this is the 
usual mean width from convex geometry if A = Aiso) as in the standard STIT setting: 

[Canonical split intensities] We have 



One further requirement reducing unnecessary generality is to drop the dependency of $s 
on s, which is referred to as time homogeneity with no need of further comments below. 
Finally, we want the split kernels $(-|-) to be isometry invariant and spatially homogeneous, 
where by the latter we mean: 

[Spatial homogeneity] We have 



We shall say that a family of split kernels with canonical split intensities, which are isom- 
etry invariant and spatially homogeneous is shape-driven and we shall stick to this level 
of generality in our considerations below. Note that the concept shape-driven as we in- 
troduced above is meant to refiect the fact that the split mechanism considered here only 




(10) 
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takes into account the splitting cell's geometry, while setting apart its size, as opposed 
to the work done so far in the literature (cf. [2]), where the focus was mainly put on 
size and split intensities whereas the shape-related aspects of splits were often reduced to 
simple geometry-independent apportionments, see ibidem and Subsection 3.3 below. A 
nested random tessellation constructed by a family of shape-driven split kernls having ad- 
ditionally the spatial Markov property is what we call a shape-driven nested Mcirkov 
tessellation in this paper. 

3.2 Kernel Regularity and Main Results 

To proceed, consider the natural space JC of compact convex bodies in R*^ and let /Co C /C be 
the subspace thereof, containing bodies centered at the origin 0, /Cq := {K G JC, c{K) = 0}, 
where the center of a convex body K stands here and below for some henceforth fixed shift- 
covariant selector of K, say the Steiner point, the circumcenter, barycenter etc. Finally, 
we let /Co;i := {K e /Co, A([i^]) = 1} be the family of centered convex bodies in M'' 
with unit split intensity, see (10), also recall from Crofton's formula [18, Thm 5.1.1.] that 
Aiso([-f^]) = 1iVi{K) with Vi standing for the order one intrinsic volume and with the 
integral-geometric constant 71 given by 

^ 2/t,_i ^ r{d/2) 

dKa ^r((d + l)/2)' 

where Kj stands for the volume of the j-dimensional unit ball. The split kernel $ induces 
a natural Markovian (i.e. probabihstic) transition kernel $® (associated forward kernel, 
shrink kernel) from /C into itself, acting by splitting a given K E )C into two cells accord- 
ing to the normalized kernel $(-|-)/|$(-|-)| and then outputting one of these cells chosen 
equiprobably at random. It is also convenient to consider the obvious modification 
(re-normahzed shifted associated forward kernel, re-normahzed shifted shrink kernel) of $® 
acting from /Co;i into itself and obtained by suitably re-sizing and shifting the outcome of 
in the unique way to make it remain in /Co;i. The obvious intermediate option, shifted 
shrink kernel <l>^ involving shifting but not re-scaling, shall also be considered in the sequel. 
The probabilistic kernels $^,$0 and <l>^^ generate natural shrink, shifted shrink and re- 
normalized shrink Markov processes on /C, /Co and /Cq i, respectively, and we denote them 

by 

Kf, K^, Kf^\ °, ... and Kf'''\ . . . 

with the law of Kf_^-^ given by ^'^{■\Kf), with Kf''^ arising as Kf shifted to the origin 
and, eventually, with Kf'^'^ arising as Kf re-sized and shifted to stay in /Co;i- Note that a 
shrink process should be regarded as a realization in law of a randomly picked branch/chain 
in the successive forward cell division procedure of Proposition 1 for the split kernel $. 
Below it will be also convenient to consider the corresponding time-reversed growth process 
modeling the backwards history restoration along a random split branch, from tessellation 
cells back to their parents in successive sub-divisions. 

A continuous-time version of the shrink dynamics will also be in focus of our interest. 
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namely we shall let Kt, t >0 be the /Cq- valued process unfolding according to the following 
stationary Markovian dynamics in time t: 

[Continuous shrink dynamics] 

[CSDl] During the time interval [t, t+dt] linearly re-size the body Kf^'^ with scaling 
factor 1 + dt, by keeping it in addition centered at to stay within )Cq. 

[CSD2] With intensity A{[Kf'^])dt do split and replace K^'^ according to the kernel 

The so-defined dynamics admits a number of stationary regimes which are, nevertheless, 
unique up to scaling under rather general regularity conditions, as we will see in the sequel. 
Proceeding toward the formulation of the core of our theory, throughout this paper we shall 
impose some regularity conditions on $ and We have decided to make them possibly 
somewhat stronger than needed, in order to keep the presentation and the argument below 
considerably simpler and more elegant without reducing the class of interesting examples 
and applications. 

[Split kernel regularity] 

[SKRl] Let for any convex polytope c C M*^, fc be a continuous density on [c] with 
respect to some hyperplane measure A <^ Ajgo satisfying fadctdH) = fc{dH) for 
any a > and /^(c) = fc for any isometry : M"^ — )■ M*^. We assume that the 
splitting hyperplane is chosen from the normalized hyperplane measure having 
this density, this is to say, the split kernels are of the form 

^dH\c) = MH)A{dH). 

[SKR2] The probabilistic transition kernel admits the unique invariant measure 
ro$ on /Co;i for the stationary re-normalized shrink process {Kf''^'^)i>i. 

Remark 1 We expect that [SKRl] already implies [SKRl], but we were not able to give 
a rigorous proof of this conjecture. 

Recall that we say that a nested tessellation {Y (s)) se[o,t] is compatible with a split kernel $ 
iff the statement of Proposition 1 holds with $5 = $ there. The main result of this paper 
is 

Theorem 1 [Existence and uniqueness theorem] Assume that ^ is a Markovian shape- 
driven split kernel satisfying the split kernel regularity conditions from above. Then for 
each t > there exists a unique nested tessellation — {Y^{s))s^[o^t] I^'^ compatible 
with The tessellation Y^ enjoys the following properties: 

(a) y*(t) is spatially translation invariant (i.e. stationary). 

(h) The tessellation satisfies the scaling relation Y^[t) = aY'^{at), a > 0. 
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(c) We have Y^{s) — Ls{Y^{t)) for t > s, that is to say is consistent in time. 

(d) The typical ce// TCell(F*(t)) ofY^it) with distribution Q*(t) has finite mean volume 
and satisfies the scaling relation 

Q*(i)-r^]0Q*, (12) 

where [a]© stands for the usual spatial scaling operation on the argument distribution 
with a positive factor a and where Q* = Q*(l). Moreover, we have 

EA([TCell(r*(i))]) = - (13) 

t 

and, letting Q*(t) stand for the image of the typical cell law Q*(t) under the standard 
re-normalizing map JCq 3 c^^ c= G /Co,i, 

Q*(l) = = ru^. (14) 



Furthermore, the integral representation 



'{t)= / {[t-'A{[c])]Qzu^)Q''{dc). (15) 



holds. 



(e) The typical cell distribution at the time 1, Q* = Q*(l), is the unique stationary law 
for the [Continuous Shrink Dynamics] satisfying (13) with t = l there. 

(f) The the law iJ^ of the zero cell ofY'^it) arises as the area-weighted modification of 
that of the typical cell, that is to say for c 3 0, 

4 Vol,(c)Q*(dc) 



As already mentioned before, the crucial feature of Theorem 1 is that the unique whole- 
space tessellation Y^{t) gets constructed in a global construction greatly reminiscent of the 
(by now) well known Mecke-Nagel- Weiss (MNW) global construction developed in [11] for 
the STIT processes. Hence, the connection between the incremental MNW-construction, 
represented here by the [Recursive split dynamics], and the corresponding g'/o&a/ MNW- 
construction, represented here by the global construction in our Subsection 5.1, carries over 
to the general class of shape-driven tessellations. It is also important to note that in view of 
our assumption on the stationary regimes for the continuous shrink dynamics, the point (e) 
of Theorem 1 combined with (13) allows to determine the typical cell Q*(i) and Markov- 
chain Monte-Carlo simulation can be an effective option. 
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Reiiicirk 2 It should be emphasized that, whereas our main Theorem 1 does guarantee 

the existence and uniqueness of the (^-compatible whole-space tessellation, it should not be 
mistakenly believed to give any kind of spatial consistency as in the STIT case. In general 
we should not expect that Y^{t) fl W coincide in law with the finite volume process we 
would obtain if we carried out the recursive split construction restricted to a proper subset 
W <Z from the very beginning. Thus, unlike many other natural properties showing 
up in our general Theorem 1, the consistency seems restricted for STIT and other closely 
related processes. 

3.3 Examples 

We shall present now some classes of examples to which our theory may be applied: 

(1) We may consider the smooth constant density a-l[H E [c]] with a > 0. This leads 
in view of (7) to a time-changed STIT tessellation y*(t). More precisely, we have by 
(2) the relation F*(t) = ^lA(t), where Yi^(t) is the standard STIT tessellation with 
driving hyperplane measure A. 

(2) For a convex set K d and r > we denote by ero(ii', r) the (possibly empty) 
erosion of K and a ball with radius r, this is the set of all interior points of K that 
have distance at least r to the boundary of K. The natural density G ero(c, r)] 
is unfortunately not smooth on s. However, for any £ > there exists by Urysohn"s 
classical lemma a continuous function /c with /c = 1 on [ero(c, r)] and /c = on 
[c \ ero(c, r — e)]. Using this density, the split kernel 

^{dH\c) := UH)HdH) 

is an admissible split-kernel in the sense of our theory, as long as A -C Ajgo. The 
restriction to the sets ero(c, r) implies that the splitting hyperplanes have at least 
distance r — £ to the boundary of their mother cell, which is some kind of hard-core 
condition added to the construction of STIT tessellations. 

An interesting class of examples can be extracted from [2], although Cowan's approach is 
entirely different from ours, as he works always in finite volumes and with a fixed number 
of cells. Ibidem, a random division rule of a cell c was called a geometry-independent 
apportionment if it only depends on the volume Vold(c) of c. Whenever a cell c splits it 
creates sub-cells c'^ and c~ and the volume proportion U — '^y^^'^c) distribution function 
G{u\v) with V — Vol(i(c). Assume that G is symmetric, i.e. G{u\v) — 1 — lim G{x\v), 

which is to say that U and ^ylfff^c) identically distributed. For normalization assume 
further G{0\v) — and lim G{u\v) — 1. We consider now a cell splitting mechanism 

with the property that U as defined above follows the distribution function G and that 
the dividing hyperplane is chosen with respect to a hyperplane measure A ^ ^iso- This 
corresponds to a geometry-independent apportionment in the sense of [2] and Theorem 1 
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ensures the existence of a whole-space nested Markov tessellation following that sphtting 
rule. Two particular examples are 

(3) G{u\v) — u, the uniform distribution on [0, 1], 



(4) G{u\v) = 1[0 < ii < ll^pVrj / f'-^l-ty-^dt, a symmetric Beta-distribution with 



^ I"; Jo 

parameter a > 0. Note that this example shares some features with Poisson-Dirichlet 
partitions considered in [1, Chap. 2.2] in the framework of fragmentation theory. 

A combination of these and other splitting mechanisms lead in view of the existence re- 
sult from above to a wide range of constrained nested Markov tessellations with rather 
homogeneous cell sizes as needed for applications to real world problems. 



4 Geometry of Shape-Driven Nested Markov Tessel- 



After we have ensured the existence of whole-space shape-driven nested Markov tessella- 
tions by Theorem 1, we turn now to some aspects of their geometry. One aim is to derive 
some general mean values and mean value formulas for this class of random tessellations 
and to make connection and to add a broad class of examples to the recently developed 
theory from [26]. The reason, why we mostly restrict ourself to the particular interesting 
cases d = 2 and d = 3 is, firstly, that with increasing space dimension, the complex- 
ity of geometric objects determined by the tessellation grows rapidly. Secondly, many of 
the mean values depend on parameters, which are determined by the precise cell splitting 
mechanism and the number of these parameters and their meaning becomes quite complex. 

Let Y = be a stationary shape-driven nested Markov tessellation in M'^ whose split 

kernel satisfies the split kernel regularity conditions and denote by X a generic class of 
convex polytopes determined by Y, for example the class V of vertices or the class C of 
its cells. 

Let Xi and X2 be any two classes of convex polytopes. We say that Xi G Xi is adjacent 
to X2 G X2 if either Xi C X2 or X2 C xi. For Xi G Xi, the number of objects of class X2 
adjacent to Xi is denoted by niXiixi). Moreover, the mean value 1^X1X2 is defined as 



where Pxi is the Palm distribution with respect to Xi, see [18] for exact definitions. 
We can interpret 1^X1X2 mean number of objects of class X2 adjacent to the typical 

object from class Xi. Furthermore, we let Ax be the intensity of objects of class X, 
i.e. the mean number of class- X-objects per unit volume and denote by J^j{X) the mean 
number of j-faces of the typical object from class X, where j ranges from to the dimension 
of the typical X-object. 




lations 
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The typical cell in the general isotropic case. Note that part (d) of Theorem 1 
imphes that EA(TCell(y*(t))) is the same for aU admissible spht kernels $ and is especially 
the same as for the typical cell of a Poisson hyperplane tessellation having intensity measure 
tA. Unfortunately, this does in general not imply that the mean intrinsic volumes of 
TCell(F*(t)) coincide with that of the associated Poisson typical cell. However, for the 
special choice A = Ajgo we have A([ir]) = ^^^^Vi{K) for any compact convex K C R'^, 
which implies 

Eyi(TCell(r*(t))) = f^j, (17) 
in particular ir/t in the planar and 6/i in the spatial case d = 3. 

The plancir case. Most of the first-order geometry of shape-driven nested Markov tes- 
sellations in can be determined by the vertex geometry of these tessellations. It follows 
directly from the construction and our assumption that the law of Y'^(t) is absolutely con- 
tinuous with respect to that of the STIT tessellation Y\^^^{t), that such tessellations can 
- in the planar case - only have T-shaped vertices, which is to say that from any vertex 
we have exactly 3 outgoing edges and two of them are coUinear. Here and in the sequel, 
we understand by an edge a line segment bounded by two vertices, but with no further 
vertices in its relative interior. In contrast, by a side of the tessellation, we understand any 
of the cell's sides. These line segments are of course bounded by vertices, but may have 
further vertices in their relative interiors. The classes of edges and sides of the tessellation 
are denoted by E and S, respectively. Moreover, we consider the class / of its I-segments 
introduced in Section 2. Observe now that any vertex of the tessellation is ... 

- ... endpoint of exactly 3 edges, 

- ... endpoint of exactly 4 sides, 

- ... endpoint of exactly 1 I-segment and 

- ... corner of exactly 2 cells. 

Using these properties, we can relate the intensities Ag, Ag, A/ and Ac to the vertex 
intensity Xy as follows: 

3 11 

A^; = 2''^^' ~ '^^v, — and Ac = 2^v- 

Moreover, denoting the mean lengths of the different typical line segments by Le, Ls and 
Lj, we get 

Li = 2Ls = 3Le, (18) 

by noting that each interior point of a side is counted twice. Moreover, application of the 
standard mean value formulas from [18] or [23] leads to /icv = 2 -|- 2 Ay /Ac = 6, i.e. the 
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typical cell has 6 vertices on its boundary in the mean. On the other hand, it is well known 
that the mean number of corners voiC) of the typical cell satisfies 



.„(C) = ..(C) = fei) = f=il=4. 

which can be interpreted by saying that the typical cell is a quadrangle in the mean with 
two additional vertices on its boundary, which are no corners of it. 

In the special isotropic case A = Aiso we can combine (17) with the fact that the first 
intrinsic volume of a convex set equals half of its perimeter length, to obtain 

p = 2Eyi(TCell(F*(t))) = — 

i 

for the mean perimeter length p of the typical cell of Y^(t). Combined with uo(C) — 4 this 
leads to precise values for L/, Ls, Le and to the observation that in the planar case the 
edge length density equals the construction time t (this is no longer true for d > 2). 
Standard mean value formulas from [23] also imply exact values for the intensities A/, As, 
As, \e and Ay. Further, the mean area a = 1/Ac of the typical cell can be found. The 
results are summarized in the following table 



La 




\e 




A/ 


Le 


Ls 


L, 


p 


a 


fJ'VE 


fJ'CV 
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TT 


TV 


TT 


iF 

TT 


TT 1 
3 t 


TT 1 
2 t 






TT 1 

2 t2 


3 


6 


4 



Note, that these mean values coincide with those for stationary and isotropic STIT tes- 
sellations in the plane, see [15], which shows that the first-order geometry of isometry 
invariant planar shape-driven nested Markov tessellations is rather rigid due to their topo- 
logical structure. However, we expect that the additionally introduced flexibility has an 
important influence on the second-moment structure of these tessellations. 

The spatial case. According to the recent theory presented in [26] , a number of the most 
basic topological mean values of a stationary spatial random tessellation is determined by 
a system of 10 parameters Hve, IJ'EP, IJ'PV, IJ'CV, IJ'Ce, IJ'CP, ^, and r, where P stands 
for the class of plates of the tessellation and by a plate we mean a 2-dimcnsional convex 
polygon bounded by edges but with no vertices or edges in its relative interior. Remarkably, 
this system of parameters can considerably be reduced in our case. To do so, we observe at 
first that the vertices of a shape-driven nested Markov tessellation in can only be of two 
types: T-vertices and X-vertices, see Figure 4. This follows, similarly as in the planar case, 
from the fact that the law of Y^(t) is absolutely continuous with respect to that of the 
STIT tessellation Y\.^^(t). Now, we can give an interpretation of the four different Greek 
parameters mentioned above. The parameter k is the proportion of T-vertices present in 
the tessellation, whereas ^ = 1, because all edges arc so-called 7r-edges in the sense of 
[26]. The parameter ■0 is the expected number of cell-side interiors adjacent to the typical 
vertex and r is the expected number of plate-side interiors adjacent to the typical vertex. 
From the geometry of the vertices we can see that any T- vertex is adjacent to exactly one 
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Figure 1: The two possible types of vertices: T-vertices (left) and X-vertices (right) 

cell-side interior and to exactly one plate-side interior. Moreover, any X- vertex is adjacent 
to exactly 4 cell-side interiors and to exactly 2 plate-side interiors, so that we get 

ip = 1 ■ it + A- {1 - k) = 4-3k, t^1-k + 3-{1-k,) = 1 + k. 

Moreover, we obviously have fiyE = 4 and fiEP = 3. Also the parameters /icy, fJ'CE and 
licp can be expressed by n and fj^py by using the standard mean value formulas for spatial 
random tessellations from [23]. Thus, the only remaining parameters in our case are k and 
X — I^PV- 

Denote by Xvit] and \v[x] the intensity of T- and X-vertices, respectively, i.e. Xv[t] = i^^v 
and Xv[x] = {1 ~ k,)Xv- We use now the vertex geometry to obtain the intensities A/, 
Xs[c], Xs[p] and of I-segments, 5'[C]-segments, 5'[P]-segments and edges, respectively. 
Here, by an edge we mean again a line segment bounded by vertices with no vertices in its 
relative interior. A 5'[C]-segment is the (1-dimensional) side of a cell and a -S'fPj-segment 
is the side of a plate (thus, the classes S[C] and S[P] are multi-sets of line segments with 
multiplicities 1 and 2 and 1, 2 and 3, respectively). Moreover, the I-segments are the sides 
of the maximal polygons inserted during the construction of the tessellation. We observe 
now that each T- vertex (X- vertex) is the endpoint of ... 

... exactly 2 (0) 1-segments, 

... exactly 4 (4) edges, 

... exactly 6 (0) -S'fCj-segments, 

... exactly 10 (8) 5* [P] -segments. 

Using these facts we get: A/ = ^ ■ 2\v[t] = /tAy, A^ = ^ ■ 4Av = 2Av, Xs[c] = ^(10Ay[T] -|- 
8Ay[x]) = (4+fi;)Ay and Xs[p] = |(6Ai/[t]+0Av[x]) = 3/tAy. Furthermore, the mean lengths 
Li, Ls[c], Ls[p] and Le can be calculated from XjLj = \eLe = Ly, Xs[c]Ls[c] = 2Lv and 
Xs[P]LsiP] — 3ivy, where Ly is the edge-length density of the tessellation, so that we end 
up with the general relationship 

Li - -Ls[C] - ~^^S[P] - -Le, 
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which in particular shows that the mean number of vertices in the relative interior of the 
typical I-segment equals whereas 

Ls[c\ , 4 - 3k Ls[p] 2- k 

/^relint(5[C])V = i- — — ^ and /Xrelint(5[P])y = i- — — T, 

Lie OK Lie K + 4 

where relint(X) stands for the relative interior of the typical object of class X. 
The geometry of the classes of objects of dimension 2 is considerably more involved and 
beside k also our second parameter x enters the expressions. Prom the general formulas in 
[26, Table 2] we obtain at first Xc — ^^Xy and with the Euler-type relation Xy — Xe + 
Ap — Ac = it follows Ap = ^Ay. Moreover, we have 

Ap = [~^^'^~ ^ J 

Note, that by a facet we mean each facet of a cell, which implies that the class F of all 
tessellation facets is a multi-set with multiplicities 1 and 2. This together with XpAp = Sy 
and ApAp = 2Sv implies now 

Ap = ; rAp. 

12 + x{k-2) 

Moreover, we can calculate t'o(-P) = ^i(-P) ^-nd ^o{F) — i'i{F), the number of corners of 
the typical plate and that of the the typical facet: 



6^^ ' "^^^^ 12 + x(«-2)' 

see [26, Table 6]. Observe now that any edge of the tessellation is contained in the relative 
interior of precisely one cell facet, whence 

_ Ap _ 2x 

/^relint(F)S " A^ " 12 + - 2) " 

Next, we note that any T- vertex is contained in the relative interior of exactly one facet 
and that there is no facet in whose relative interior an X-vertex is contained, which implies 

- l2 + x(«:-2) 

and the general relationship K//rciint(F)£; = 2yUi.ciint(F)y- Moreover, jipp = 2Ap/Ap = 
12/(12 + — 2)), because each plate is part of two facets. Unfortunately and in contrast 
to the plates and facets, the mean values for maximal polygons cannot be expressed by the 
two parameters k and x and we skip their discussion for this reason. For the approach in 
the STIT case we cite [24] and remark that the method developed there cannot be extended 
to any other shape-driven nested Markov tessellation. 
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To derive formulas for the typical cell C of the tessellation y, we can utilize the standard 
mean value formulas for spatial tessellations from [23, Chap. 10.4], which lead to 

fJ'CV — 7^ ! ^J'CE — 7^ ! fJ'CP 



6-x' 6-x' 6-x 

and furthermore the formulas from [26, Table 6] to obtain 

6-x 6-x 6-x 

for any shape-driven nested Markov tessellation in M^. 

From our results, we can deduce some inequalities for the basic parameters k, and x- First, 
observe that from 

i.o(P) = |(4 + «:)>^X>3 

we get X ^ |- Moreover, we obviously have x < 6, because iicp for example has to be 
positive. Moreover, iyo{C) > 4 implies 

^ 12 - 2x 

X 

Remcirk 3 It would be interesting to know if x can be arbitrarily close to 6 and if there 
are examples for which n becomes arbitrarily close to 1. These questions are closely related 
to open problems raised in [26]. 

The S TIT- case. Mean values for planar and spatial random STIT tessellations were 
considered in [15], [16] and [24]. Using the additional iteration stability and the property 
that STIT tessellations have Poisson-typical cells, the parameters Ay and La in the planar 
case and Ay, k, x, Ly and Sy in the spatial case can be made available. In the planar 
case, the two parameters are La = t and Ay = t^C, with C, given by 



C = / / [ui,U2\ll{du2)Tl{dui), 

JSi JSi 



where IZ is the directional distribution introduced in Section 2 and where [mi,M2] stands 
for the area of the parallelogram spanned by ui and U2. In the spatial case, we derive from 
TZ the two constants 



C2 = / / [ui,U2]Ti{du2)Ti{dui), 
J S2 J S2 

C3 ^ [ui,U2,U3]TZ{du3)7l{du2)TZ{dui), 

J 09 S2 "J 09 



'1S2 ''52 >'52 

where, similarly to the planar case, [-Ui, ""s] is the volume of the parallelepiped spanned 
by Ml, U2 and 7/3. By a direct verification it can be shown that in the isotropic case, C2 and 
(■3 are given by | and | , respectively. We are now in the position to express the parameters 
discussed above by the construction time t and the two numbers ^2 and 



t, Ly — 5'y(2, Ay — 5'yC3, — X 



2 36 
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5 The Global Construction and Proof of Theorem 1 



As already signaled above, our proof of the principal Theorem 1 goes by providing first a 
global construction for the infinite volume ^-splitting tessellation, constituting a direct and 
natural extension of the procedure given by Mecke, Nagel and Weiss in [11, 12] for STIT 
tessellations. We then show that the so-obtained tessellation is the unique one compatible 
with $ in the full space M"^ and we verify a number of properties thereof, obtained as 
by-products of the argument. 



5.1 The Global Construction 

The idea underlying our global construction involves a gradual reconstruction of the space- 
tessellating field, given its parts already reconstructed, starting from the void and eventu- 
ally covering the whole of M.^. To this end, we consider first the time-reversal of the shrink 
processes and the associated Markovian (probabilistic) growth kernels. For k,k' e JCq-i we 
let 

^*-''-'\dk\k') = FiK^-''-'' e dk\K^','' = k') = ^^''''^^lf]j^^^^^ (19) 

which is defined ro^-almost everywhere and is referred to as the re-normalized growth 
kernel in the sequel. Dropping the re-normalization we also define, with k, k' e /C, 



^*{dk\k') = mf e dk\Kr_, = k') = ^^idm-^-mM ^^o) 

zu^{d{k' )o;i) 



where (A;)o;i stands for the standardized version of k, re-sized and shifted so as to fall into 
/Co;i- Below we grant the name growth kernel to the so-defined $*. Given the Markovian 
growth kernel $* and a seed cell Sgccd G /C we consider the spinal growth process 
Sq = Sseed, Sl,S2, . . . , Corresponding in law to the time-reversed shrink process with infinite 
indexation, satisfying 

P(5^i e dk\S* = k') = ^*{dk\k'). (21) 

For the so-defined process, we clearly have by definition S*_^_i D S* almost surely for 
i = 0, 1,2,. . . . 

To make proper use of the spinal growth process we shall now ascribe auxiliary time marks 
S* T* to its subsequent cells, with the objective of having the joint process {S*,r*)i>o 
described in law as follows: 

• S* - the subsequent ancestor cells, backwards in time along a randomly picked branch 
of the recursive split procedure generating the tessellation considered, 

• T* - the corresponding cell's birth times in the course of the recursive split procedure, 
where by birth time r* of cell S* we understand the time moment in [0, t] , where the 
super-cell S*_^_i of S* splits and gives rise to S*. Clearly, < . . .t*^i < r* . . . < t 
almost surely. 
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By convention, we put Tq — t for the seed cell. Then, aiming at the above interpretation 
and recalling (10) determining the total split intensity = ^([/cj]) for a cell ki, 

we conclude that we have, conditionally on S*, S*_^_i and r* given, 

because the corresponding conditional probability of surviving the time period [ti+i,ti) in 
state ki upon shrinking to it at tj+i is just A([A;i])e~^*^'*^*])(**~*'+i)dti+i by the construction 
introduced in Section 2. Solving the proportion yields 

P(t*,i e dti+,\S*_,, = S* = r* = ti) = A([fc,])e-^(['='])(*'-*^+i)di,+i. (22) 

Consequently, by (22) above, the extended Markovian kernel 

^*{dk,dt\k',t') = <!>*{dk\k')A{[k'])exp{-A{[k']){t' -t))dt (23) 

governs the joint Markovian dynamics of the time-marked spinal growth process {S*, T*)i>o. 
We close this paragraph by noting that the cells 5"* eventually cover the whole space, which 
is Ui^o ~ "^^^^ '^^^ ^^^^ follows, see also [11, Lem. 4.1]: Denote by a < t the 
largest birth time of a cell from the spinal growth process hitting the convex hull of Bi (p) 
and the center c{Ssccd) of the seed cell, where Bi{p) is the ball around p G M"* with radius 1. 
Then p is contained in the interior of S** for any r* < a, thus the event {3i : Bi{p) C S**} 
has probability one. The claim follows now from the fact that M.^ can be covered by balls 
Bi{p) form a countable set of midpoints p and the fact that the process S* is growing, that 
is S* C S*^,. 

The Construction. The following crucial procedure is our global construction with 
initial seed cell Sgeed, outputting the tessellation Yg^^^^{t), t > 0. It can be regarded as a 
direct generahzation of the global construction of STIT tessellations from [11, 12]. 

[Construction: Spinal phase] Starting from the initial seed cell S'seed at terminal time 
t, construct the time-marked spinal growth process {S*,T*)i>o backwards in time as 
determined by the extended spinal growth kernel $*. sec (23). In this way a (non- 
stationary) frame tessellation of the entire is obtained, consisting of an infinite 
chain of nested cells fully covering the space. 

[Construction: Local fill phase] Given the spinal chain with time-marked cells, for 
each cell c, having time mark falling below the threshold t, carry out within c the 
usual [Recursive split dynamics] with split kernel $ up to time t. Output the 
resulting tessellation as Ys^^^^ (t) ■ 

5.2 Completion of the Proof 

The further argument in our proof of Theorem 1 relies on showing that the desired unique 
tessellation y* can be obtained from the global construction from the last subsection 
under appropriate (randomized) choice of the initial seed cell Sseed made coincide with the 
corresponding zero cell for F*. The reasoning splits into several steps and goes as follows. 
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The existence of a translation inveiriant field. We identify a random tessellation 
in M*^ with the union of its cell boundaries, i.e. we regard it as a random element in the 
space T of closed subsets of M^. For our split kernel satisfying our imposed regularity 
assumptions, we consider the family 

J^={£(y*(i,[-i?,i?]<^)): i?eN} 

of laws of tessellations [— -R, i?]), G N. As the space T is compact, the familiy y is 

tight. Whence, by Prohorov's theorem, there exists a convergent subsequence in y ^ whose 
limit is the distribution of a random closed set Y^{t\ oo) in M*^. We note that Y^{t\ oo) is 
the frame of a random tessellation, the cell interiors of which are the connected components 
of l^*(t; oo)'s complement. Moreover, go) is ^-compatible, as obviously the kernel $ 

has the Feller property due to our assumption [SKRl] on the densities /g. Consider next 
the family 




of versions of Y{t\ oo) smeared (shifted and normalized) over increasing cubes. Using again 
the tightness and passing to a subsequence we find in the limit F*(t), which is ^-compatible 
and translation invariant. Moreover, we readily conclude from the construction that 
can be taken consistent in time and that the desired standard scaling relation 

Q;y*(Q;i)=y*(i), a>0, (24) 
holds in view of (10), possibly upon further passing to a subsequence. 



The typical cell and the zero cell. Write A*(t) for the cell density of Y^{t) and 
let Q*(i) be the corresponding typical cell distribution on /Cq, so that the standard cell 
intensity measure decomposition [18, Eq. (4.2)] becomes 

e*(i) = AS (t) Shift (£''®Q*(t)) 

where 0*(t) = ^^c^Q^w^iy^ity^^c is the particle intensity measure in the sense of [18, 
Section 4.1], i'^ is the usual d- dimensional Lebesgue measure, whereas Shift(-) stands for 
transporting the argument measure via the shift map [x^c] ^ x ^ c with a; e M'^ and 
c e /Cq. Under this notation, the [Recursive split dynamics] of Proposition 1 yields in 
view of (10) 

j\%{t) = \%{t) ^ kiXcWmc) = AS(t)EA([TCell(y*(t))]). 

Indeed, this follows readily from the definition of the generator together with the definition 
of the typical cell, see [18, Eq. (4.6) and Thm. 4.1.3]. However, by stationarity and by 
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the scaling relation (24), there exists a constant a > not depending on t such that 
A*(i) = af^. Inserting this expression into the last equation immediately leads to 
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A([c])Q*(t)(dc) = EA([TCell(F*(t))]) = (25) 

/Co 

which was part (13) of statement (d) in Theorem 1. In particular, this shows that the 
size of the typical cell of Y'^{t) - measured here by A([-]) - is influenced only by the 
split intensities (4), i.e. the total masses of the split kernel $, and not by the precise 
split geometry determined by it. Moreover, it ensures in view of the well-known Usryson 
inequality from convex geometry (cf. [8, Thm. 8.9]) that the typical cell has finite mean 
volume, adding to the statements of part (d) of the Theorem. 
By the same argument as used above, again by (8) and (10), we have 

^^-e^it) - AS(^)Shiftf£'^^ / [2^^{-\c)-SM[c])Q''mdc] 



dt ' ' ' V 



/Co 



= AS(t)Shift ® E[2$0(.| TCell(y*(t))) - 5TCeii(y*w)]) , 

where the inner bracket reflects the operation of replacing the splitting cell c by the two 
new cells arising. Putting the equations for A*(t) and ©*(t) together yields 



±0^ 
at 



it) = I [2$^(-|c) -5,- q'{t)]K{[cW{t){dc) (26) 



for the typical cell distribution Q*(t). Observe next that the use of (24) gives us the 
important scaling relation 

Q*(i) = r^]0Q*(l), (27) 
whence the desired equation (12) follows upon putting Q* := Q*(l). 

To proceed, we let now Q*(t) stand for the image of the typical cell's law Q*(i) under the 
standard re-normalizing map 1Cq3 a-^ c= -j^^^c G ATq,!. Hence 

(f{t) = Q* = const. (28) 

Consequently, (26) leads us to 



0= / (2^^Mc)-5c-(f)d(f{dc) 



and thus 

$o^,(-|c)Q'^(dc), 

/Co,i 

SO that eventually, but not unexpectedly, turns out to coincide with the unique invariant 
measure zu^ for the re-normalized kernel ^q i, i.e. 

= ro$ (29) 
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which is claimed as (14) in (d). 

Getting back to the original non-rescaled typical cell Q* constituting the focus of our 
present interest, from (29) we obtain the integral representation 

Q*(t) = / {s(Dw<s,)u{ds\t), (30) 

where u{ds\t) = Q*(t)({A([c]) G ds}) = P({A([TCeU(F*(t))]) G ds}) is the law of the 
typical cell's split intensity /mean width under Q*(i). However, (15) is just a reformulation 
of (30). 

The corresponding statement from part (f) for the zero cell of the tessellation follows now 
readily by the general theory, see e.g. [18, Thm 10.4.1], because of the finite mean volume 
of the typical cell, which was part of statement (d) of the main Theorem. 
Further, in part (e) the claimed preservation of Q* by the [Continuous Shrink Dynam- 
ics] follows now by 

• using (27) to write [1 + dt] Q* = Q*(l - dt) and comparing with [CSDl], 

• noting thereupon that during the time interval [1 — dt, 1] of the incremental MNW- 
construction, when constructing y*(l) out from y*(l — dt), precisely the contents 
of rule [CSD2] is carried out for each cell of y*(l — dt). 

Alternatively one can just use (26) at t = 1. Finally, the required uniqueness of the sta- 
tionary regime of the [Continuous Shrink Dynamics] follows from (13) together with 
the observation that under each of the conditions EA([i^f '°]) = a > there is precisely 
one stationary regime by assumption [SKR2]. 

Identifying the translation invciriant field as the outcome of the global con- 
struction. To proceed with our argument, consider the sequence {S*,T*)i>Q, S*^i D S*, 
of successive time-marked cells of y*(t) containing the origin 0. Then it is clear that 
{S*, T*)i>o is a realization of time-marked spinal growth process backwards in time as de- 
termined by the extended spinal growth kernel given by (23) and with initial seed cell 
Sq — Sseed given in law as the zero cell (16) of Y'^{t). Thus, Y'^{t) admits a spinal frame 
arising as an instance of [Construction: Spinal phase] initiated by the zero cell (16). 
Moreover, once this spinal frame for F*(t) is given, the required compatibility of F*(t) 
immediately imphes that the remaining part of the tessellation has to coincide with the 
outcome of an instance of our [Construction: Local fill phase] as required. 

Uniqueness in the translation invariant regime. This follows directly by the pre- 
vious paragraph where our identification of Y'^(t) with the outcome of global construction 
was only reliant on the fact that y*(t) is ^-compatible. 
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Uniqueness in the general regime. Write Zf for the zero-cell of Y^{t) and set (f to 
be its corresponding birth time in the course of the recursive split construction. Further, 
let recursively '^t+i be the parent cell of Z* and, again, its corresponding birth time. 
Next, consider some arbitrary ^-compatible nested tessellation Y'(t) and define {Z'i^,('k) 
for Y'{t) in analogy to (Zf,C|) for Clearly, both {Z'^X'k) and (Zf,^) are Markov 

processes governed by the same Markovian dynamics driven by suitable growth kernels, 
recall (23) in Subsection 5.1. In view of the assumed [Split kernel reguleirity] we conclude 
that 

hm dTVar((Z'„a),(Z*,C|'))=0, 

with dxvar standing for the total variation distance of distributions, cf. [9, Thm. 20.20]. 
Consequently, given e > the processes l^'(t) and can be coupled so that with 

probability exceeding 1 — e we have Z* = Z'^ and = CI: as soon as A; is large enough, 
say k > k[e\. We modify this coupling by erasing, on the above event of equality, both 
Y'(t) and Y^(t) over the common region Z* = Z'^, and then restoring both over this 
region using the same instance of the recursive split construction. The resulting new 
coupling has the property that, with probability not smaller than 1 — e, we get the equality 
Y'{t) n Z'fc = Y^{t) n Z'fc for k large enough. RecaUing that almost surely Z'^ t shows 
that dTVar(^'(^), ^*(^)) < e- Taking into account that e > was arbitrarily small we finally 
conclude that Y'{t) and Y^{t) coincide in law, as required. □ 

A Argument for Proposition 1 

In order to keep the paper self-contained, we sketch here the proof of Proposition 1 and 
refer to Lemma 1 of [19] for more details in a more general context. 

We start by defining for y e Y[o,t](VF) the split measure on [1^] x 1^ by 

^Iv'-^ / ^{H,McnH))^Y{dH\c,y), 

ceCells(j/) "^W 

where ttq stands for the metric projection of the origin onto the argument closed convex 
set and denotes the unit mass Dirac measure. Write Facets (x, i/) for the facet that 
arises by splitting by H at time s the unique cell c G Cells(|/) with 7io{cn H) = x if such a 
facet exists and put Facets (a;, H) $ otherwise. Then the generator can be expressed 
as 

[^YfM- I [f(yUFe.cets{x,H)) - f(y)]Tl^(dH,dx). 

J[W]xW 

Taking into account that Y{s, V) arises as restriction of Y{s, W) to V for all compact and 
convex ^ C 1^, we get 

[h^f]{y) = E[[h^f]{Y{s, W))\Y{s, W)nV^y]. 
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Let be the cr-field generated by C R'^, s e [0, t]. Defining the -measurable 

random measure T^w by T^w :— T,^' , we conclude from our discussion above that for 
any measurable B G [V] x V we have 



ensures the existence of an -measurable random measure on x M satisfying 



with = T^^^*^ and we readily conclude that the generator of the (non-homogeneous) 
Markov process Y{s), s e [0,t] is given by 



where the last expression is defined for all bounded /. which are -measurable for some 
bounded convex W C M.'^. The remaining relation (9) is an easy consequence of our above 
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